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Probabilities for two-neutron transfer reactions, P 2 n, are often discussed in comparison with the 
square of the corresponding probabilities for one-neutron transfer process, (Pin)^, implicitly assum¬ 
ing that (-Pin)^ provides the probability of two-neutron transfer reactions in the absence of the 
pairing correlation. We use a schematic coupled-channels model, in which the transfers are treated 
as effective inelastic channels, and demonstrate that this model leads to P^n ~ (Pin)^/4, rather 
than P 2 n = (Pin)^, in the pure sequential limit. We argue that a simple model with spin-up and 
spin-down neutrons in a single-particle orbital also leads to the same conclusion. 
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It has been well known that the pairing correlation 
enhances cross sections for the two-neutron transfer pro¬ 
cess as compared to those in the uncorrelated limit [ll-Q ■ 
Those cross sections are often converted to the transfer 
probabilities by dividing them by the Rutherford cross 
sections, and are plotted as a function of the distance of 
the closest approach, D, for the classical Rutherford tra¬ 
jectory. This representation in fact provides a convenient 
way to discuss the reaction dynamics, since the cross 
sections for different values of incident energies and the 
scattering angles can be analysed in a unified way. The 
enhancement of the two-neutron transfer process has cus¬ 
tomary been discussed by taking the ratio between P2n 
and (Pin)^ SM , where Pin and P 2 n are the proba¬ 
bilities for the one- and two-neutron transfer processes, 
respectively. That is, it has been usually believed that 
the quantity (Pin)^ provides a reference probability for 
the two-neutron transfer process which would be realized 
in the absence of the pairing correlation [l^ . 

In this paper, we discuss the validity of this as¬ 
sumption. To this end, we consider the two-neutron 
transfer probability in the no-correlation limit, where 
the two-neutron transfer process takes place in a com¬ 
pletely sequential manner. This work is partly moti¬ 
vated by a recent result of a time-dependent Hartree- 
Fock (TDHF) -|- BCS calculation, which shows that the 
ratio P2n/(Pin)^ in the absence of the pairing correlation 
is well parametrized as [l^ . 

P2n N^-1 Nf-1 
(Pi„)2 ^ 2N, ' Nf ^ ^ ’ 

where Ny and Nf are the number of valence nucleons 
and the number of available states in the receiver nu¬ 
cleus, respectively. This equation suggests that the ra¬ 
tio P2n/(^in)^ is not Unity in general, but is more com¬ 
plex and never exceeds 1/2. In this paper, we employ 
the coupled-channels approach to investigate this prob¬ 
lem from a different perspective. In particular, we use a 
schematic coupled-channels model for two-neutron trans¬ 
fer, and attempt to understand the result of TDHF. 


In the coupled-channels approach to transfer reactions, 
one often treats transfer channels as effective inelastic ex¬ 
citations [13-123 ■ i-hi® paper, we use the same treat¬ 
ment for the transfer channels and consider the following 
coupling matrix for a sequential two-neutron transfer re¬ 
action [23, Hill 

/ 0 F{r) 0 \ 

1/ = F{r) -Q F{r) . (2) 

\ 0 F{r) -2Q j 

Here, we have assumed that all the channels have zero 
angular momentum. In this equation, F{r) is the form 
factor for the coupling between the entrance (On) and 
the one-neutron (In) transfer channels, while Q is the 
Q-value for the In-transfer reaction. In this coupling 
scheme, the On channel is coupled to the In channel, 
which is sequentially coupled to the two-neutron (2n) 
channel. The no-correlation limit is simulated by setting 
the coupling between the In and the 2n transfer channels 
to be the same as that between the On and the In trans¬ 
fer channels, and also by setting the Q-value for the 2n 
transfer channel to be exactly twice the Q-value for the 
In channel. The direct coupling between the On and the 
2n channels is also set to be zero. 

We apply this model to the ^°Ca-|-®®Zr reaction, for 
which the experimental transfer cross sections have been 
reported in Ref. [l2j| . To this end, we use a function 
which asymptotically has an exponential form, 

F{r) ~ (3) 

a 

for the coupling form factor F{r), and set the transfer 
Q value to be Q=0 [13 US]' ^^e actual calculations 
shown below, for a numerical reason, we use a derivative 
form of the Fermi function with the parameters /?, R, and 
a.) With the Woods-Saxon type for the nuclear potential, 
with the parameters of Vq = 140 MeV, tq = 1.1 fm, and 
oq = 0.65 fm for the real part and Wq = 30 MeV, rw = 
1.15 fm and aw = 0.1 fm for the imaginary part, we vary 
the parameters in the coupling form factor, Eq. m , so 
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FIG. 1: (Color online) The transfer probabilities, defined as 
the ratio of the transfer cross sections to the Rutherford cross 
sections, for the + ®®Zr reaction. These probabilities 

are plotted as a function of the distance of the closest ap¬ 
proach, D, of the classical Rutherford trajectory. The dotted 
and the solid lines denote the one- and two-neutron transfer 
probabilities, respectively, while the dashed line is a quarter 
of the square of the one-neutron transfer probability. The 
experimental data are taken from Ref. [l3| . 


that the experimental data for the one-neutron transfer 
reaction can be reproduced. To this end, the coupled- 
channels equations are solved using a version of the com¬ 
puter code CCFULL [23. The resultant values for the pa¬ 
rameters are /3 = 9 MeV fm, R = 1.15 x(40^/^ -I- 96^/^) 
fm, and a = 1.3 fm. 

Figure 1 shows the transfer probabilities so obtained. 
Here, the transfer probabilities are defined as the ratio 
of the transfer cross sections to the Rutherford cross sec¬ 
tions, that is, Pxn = (doxn!d^)/[du-n/<Kl), where x = 
1, 2 is the number of transferred neutron, daxn/dO, and 
dan/dO, are the transfer and the Rutherford cross sec¬ 
tions, respectively. This definition is applied both to the 
experimental data and to the theoretical calculations. 
The transfer probabilities are plotted as a function of 
the distance of the closest approach, D, of the Ruther¬ 
ford trajectory for the scattering angle of 0c.m. = 140 
degrees in the center of mass frame. The dotted and the 
solid lines denote the transfer probabilities for the In and 
the 2n channels, respectively. While the In probabilities 
are well reproduced, as expected, the 2n probabilities 
are largely underestimated by this calculation. One can 
clearly see that the 2n probability, P 2 n, is consistent with 
a quarter of the square of the In probability, (Pin)^/4, 
which is denoted by the dashed line in the figure. 

This relation can be easily understood if one uses 
the time-dependent perturbation theory. In the semi- 
classical coupled-channels approach, one assumes a clas¬ 
sical trajectory r(t) for the relative motion between the 
colliding nuclei, and solve the time-dependent coupled- 
channels equations for the intrinsic motion @,[ 3 . Ap¬ 
plying the first and the second order perturbation theory, 
the amplitudes for the one- and the two-neutron trans¬ 
fer processes for the sequential two-neutron transfer cou- 
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FIG. 2: (Golor online) A schematic model for the two-neutron 
transfer process of spin-up and spin-down neutrons. P and T 
denote the projectile and the target nuclei, respectively. 


pling, Eq. (g]), read. 


“ aj 

pOO 

1 dte-^^^/^F{r{t)), 

— 00 

(4) 

= (s 

\ 2 poo 

/ dte-^^*/^F{r{t)) 

/ J — 00 



dt' /^F{r{t')), 

— 00 

(5) 

1 

2 

1 "1 ^ 

- / dte-^^^/^F{r{t)) , 

J — 00 

(6) 


respectively. The last equality is due to the property of 
the pure sequential transfer, that is, Q 2 n = ‘^Qin and 
P(ln-2n)=F(0n-ln). By squaring these equations, one 
obtains P 2 nl{Pinf = |a 2 nP/|ain|^ = 1/4, which is in¬ 
deed realized in Fig. 1 for large values of D, at which the 
perturbative treatment is justified. 

The factor of 1 /4 can also be obtained with a more mi¬ 
croscopic model, which is illustrated in Fig. 2. Here we 
consider a transfer of spin-up and spin-down neutrons, 
which initially occupy a single-particle state in a projec¬ 
tile nucleus (see the state (a) in Fig. 2). One of those 
neutrons is initially transferred to a target nucleus (the 
states (b) or (c), depending on the spin of the trans¬ 
ferred neutron), which is followed by a transfer of the 
other neutron to the target nucleus (the state (d)). We 
assume that the matrix elements for the transfer process 
do not depend on the spin of the transferred neutron, 
and that the spin flip does not occur during the transfer. 
We thus have (fp \V\ fr) = (ip \V\ It) = F{r) and 
(tp \y\ ir) = (ip \V\ tr) = 0, where V is the opera¬ 
tor which induces the transfer, and P and T denote the 
projectile and the target nuclei, respectively. We again 
use the time-dependent perturbation theory in order to 
evaluate the transfer probabilities. For the one-neutron 
transfer probability, there are two distinguishable final 
states, (b) and (c) in Fig. 2, and one has to add the 
probabilities for the processes (a)—>'(b) and (a)—>'(c). One 
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thus obtains (see Eq. dU), 


Pin = 2 X 


ih 




with different values of jz, that is, the z-component of 
the single-particle angular momentum in the receiver nu¬ 
cleus. If one neglects the spin-flip components, however, 
the formula would become 


For the two-neutron transfer process, there are two indis¬ 
tinguishable paths, (a)—>'(b)—>'(d) and (a)—>-( 0 )—>'(d), to 
the final state, and one has to add the amplitudes first. 
This leads to (see Eq. dH])), 


^2n — 


1 1 ^ 

— dte-^^*/^F(r(t)) x2 


( 8 ) 


Comparing Eq. 0 with Eq. dSl), one again obtains 
P 2 n/{Pinf = 1/4. It is easy to confirm that this relation 
still holds even if one considers the anti-symmetrization 
of each state, e.g., |a) = (| fpip) ~ I 4 'Ptp))/’\/2 and 

1^) = (I tpir) - I iTtp))/V^- 

As in the multi-phonon couplings in the coupled- 
channels approach (^ . , one can make a relation be¬ 

tween the coupled-channels model of Eq. m and the 
schematic model of Fig. 2. That is, by introducing a 
single effective one-neutron transfer channel defined by 
|ln) = (|6) + |c))/y2, it is easy to find (ln|E|a) = 
{d\V\\n) = where |a), |6), |c), and \d) are the states 
shown in Fig. 2. Therefore, identifying F = \/2F, the 
two models are actually equivalent to each other. Notice 
that the other combination of the states \b) and |c), that 
is, (16) — |c))/-\/2, couples neither to |a) nor to \d) and is 
decoupled from the model space. 

The factor of 1/4 for the relation between P 2 n and 
(Pin)^ is consistent with the previous result of TDHF, 
Eq. (Il])i if one disregards the dependence on Nf. No¬ 
tice that the Nf dependence in Eq. © was obtained by 
counting the number of possibilities to put nucleons in 
the final single-particle state [l^. To this end, the prob¬ 
ability was assumed to be the same for all the final states 


P2n N, - 1 

(Pi„)2 ^ 2iV„ 


(9) 


with which one obtains P 2 n/{Pin)^ = 1/4 for Ny = 2. 

In summary, we have investigated the two-neutron 
transfer reactions in the no-correlation limit. To this 
end, we have used a schematic coupled-channels model, 
in which the transfer channels are treated as effective in¬ 
elastic excitations. We have shown that the probability 
of two-neutron transfer process, P 2 n, is given by a quar¬ 
ter of (Pin)^, that is, P 2 n/(Pin)^ = 1/4. This result is to 
some extent consistent with the result of time-dependent 
Hartree-Fock calculations for two valence neutrons. The 
two-neutron transfer probabilities have customary been 
compared with (Pin)^, rather than (Pin)^/4. Of course, 
many experimental data are for inclusive processes, and 
the enhancement factor for the two-neutron transfer pro¬ 
cess reflects not only the pairing correlation but also 
the phase space factor for the intermediate and the fi¬ 
nal states. Nevertheless, there is no strong reason why 
the two-neutron transfer probability should be compared 
with (Pin)^, and we advocate using (Pin)^/4, which has 
a clearer physical meaning as a reference probability, at 
least for a core-l-two-neutron system. 
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